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INTRODUCTION 


This report contains a compilation of the data collected by a 
combined NASA/University of Chicago field program on Siple Coast, 
Antarctica from October 1983 to January 1987 and a description of 
the methods used to obtain the data. In its first three field 
seasons, beginning in October 1983, the project has concentrated 
on studying Ice Streams B and C (see Figure. 1) . A number of 
publications have resulted from the analysis of these data and 
they are referenced in Appendix 4. This report will assist in 
the continuing analysis of these data and adds them to the 
growing data base that exists for the region. 

The NASA/University of Chicago program is part of a larger 
project which also involves the Byrd Polar Research Center at 
Ohio State University and the Geophysical and Polar Research 
Center at the University of Wisconsin. The major goals of this 
larger project are: 

1) To measure the mass balance and configuration of the major ice 
streams in the Ross Sea drainage basin of the West Antarctic ice 
sheet ; 

2) To determine the forces that control the flow of these ice 
streams ; 

3) To investigate long-term trends in the drainage basin. 

To meet these goals, the .NASA/University of Chicago group 
designed a data acquisition operation with the following specific 
objectives: 

1) Establish a network of stations to measure the surface 
velocity and strain rate profiles across the widths of Ice 
Streams B and C at locations near their grounding lines to 
determine net ice stream discharge rates. 
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2) Establish stake schemes to measure longitudinal velocity and 
strain rate profiles to determine the spatial gradients of both 
ice transport and the forces which restrain the ice stream 
motion. 

3) Establish a network of stations which surround Crary Ice Rise 
(Figure 1), to measure velocities and strain rates, and thus 
determine its influence on the flow of Ice. Stream B, and the 
state of mass balance of the ice rise. 

4) Establish a network of stations which surround a feature named 
Ice Rise "a' 1 , a feature seen on early airborne photography in the 
region, to determine its position and characteristics. (We state 
here that Ice Rise "a” is not a true ice rise at all, as its 
velocity is the same as the surrounding ice, which is itself, 
lightly grounded.) 

5) Establish a network of stations to measure the ice velocity 
and strain rate fields in the lower, lightly grounded regions of 
Ice Stream B referred to as the "ice plain". 

6) Map the precise location of the grounding lines of Ice Streams 
B and C. 

7) Perform regional studies of small-scale ice rumples to 
characterize their effect on the large-scale flow. 

8) Establish stake schemes and carry out resection surveys to 
measure the velocity field in the severely crevassed margins of 
both ice streams and ice rises, to determine the local stresses 
in this region. 

9) Develop a new method for rapid acquisition of strain rate data 
using multi-leg rosettes. 

Although the data presented within this report are mainly the 
result of the collaborative effort of the University of Chicago 
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and NASA, the wider collaboration with Ohio State University and 
the University of Wisconsin is also evident, particularly in the 
map figures. We have used the results published in Shabtaie and 
Bentley (1986) as the basis for the boundaries of the ice streams 
and their grounding lines and the radar definition of ice rises 
and rumples. We are also grateful for vigorous continuing 
discussions between the authors and the other collaborating 
institutes; in particular, with S. Shabtaie, C. R. Bentley and 
I. M. Whillans. 
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Figure 1 : Regional map showing station positions. 
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SECTION IS ICE VELOCITY 


During this program, 26 ice velocities were measured by 
accurately determining the position of stakes using a Doppler 
satellite tracking system at several times and then calculating 
the distance moved during the time intervals. To obtain the 
distance between two geographic positions, Clarke's formulae were 
used, (see Jackson 1980) , with a value for the semi-major axis of 
6378.16 km and a flattening of 1/298.25. 

A Magnavox MX1502 geoceiver was used for precise position 
determination. Each successful tracking of a satellite pass 
provides an independent fix of the receiver position. Successive 
fixes are averaged together by the receiver so that after about 
30 passes, the position reaches a steady value within a sphere of 
about a 1.6- meter radius. However, this value of the standard 
deviation cannot be considered the error in absolute position 
because there is also an error in the broadcast satellite 
ephemeris which can cause an additional error of up to 20 meters 
(Whillans, personal communication). This imprecision can be 
improved with post-processing by using either precise ephemerides 
or translocation techniques. To obtain 30 passes at 84 degrees 
south, the average latitude of the field program, the geoceiver 
must operate for about 24 hours. For operational reasons, a 
geoceiver was rarely allowed to track more than 40 passes. More 
often, only about 10 passes were used to obtain a position. 

We have adopted the error analysis of Thomas et al (1984) in 
which for less than 30 passes, they take the root-mean-square 
radial error as 17.48/Vn meters, where n is the number of passes 
used to obtain a point position. This is twice the error found 
during the geoceiver test program and about 1.5 times larger than 
the error calculated by the Magnavox MX1502 , but we feel this 
larger error should be used, given the lack of any other control. 
This error indicates that for ice moving at a velocity greater 
than 150 m/a, and a measurement interval of at least 1 year, as 



few as 10 passes will give the velocity to within 5 percent; 30 
passes will reduce the error to less than 3’ percent. For ice 
moving 100 m/a, 30 point positions are needed to give 5 percent 
accuracy, while for ice moving at 30 m/a, translocation methods 
are needed to obtain better than 10 percent accuracy. 

Table 1 summarizes the velocities measured by geoceivers. It 
includes the latitude and longitude of each station, the date at 
which the original position was observed, the number of passes 
for the two point position determinations, and the velocity 
magnitude and direction along with an estimate of their errors. 
Our computed velocities and those of the Ross Ice Shelf 
Geophysical and Glaciological Survey (RIGGS) program are plotted 
in Figure 2. 

Table 2 gives the velocity across a 54-km width of Ice Stream B, 
including station DNB. The width of the glacier at this point is 
84 + 3 km taken from the map in Shabtaie and Bentley (1987) . The 
velocities were determined by transecting a line of stakes twice 
to determine their relative positions and motion, and then adding 
the velocity determined by Doppler satellite positioning at 
station DNB. The observations were made over two 1-year time 
intervals; stakes All to B10 were surveyed in 1983 and 1984, and 
stakes Bll to B19 were surveyed in 1984 and 1985 so all the 
observations were reduced to two epochs. The rotation of the 
stake line (about 6 minutes in 2 years) was included in these 
calculations. The velocity at stake B18 was also determined by 
Doppler satellite positioning. The misclosure between the two 
methods is 5 meters in magnitude and 5 degrees in azimuth, which 
corresponds to a 26-meter error for both positions each time they 
were measured. Figure 3 plots the transverse velocity profile. 
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ORIGINAL PAGE IS 
OF POOR QUALPTY 


TABLE 1 VELOCITY MEASUREMENTS 


Site Latitude Longitude No. of Velocity Date 

Passes Magnitude Azimuth* 

deg/min/sec deg/min/sec YR1 YR2 (m/a) (deg) D/M/Y 

Crary 


E2.3A 

83 

33 

09 

.1 

168 

12 

43 

17 

9 

172 

+ 7 

319 + 

2 

30 

.1 

11 

84 

E3 

83 

22 

44 

.0 

169 

34 

33 

6 

30 

21 

8 

144 

21 

19 

. 7 

11 

84 

E4 

83 

9 

20 

.2 

171 

36 

38 

15 

7 

7 

8 

23 

49 

19 

.9 

11 

84 

K1 

83 

10 

24 

.3 

168 

9 

9 

5 

8 

225 

10 

326 

3 

20 

. 1 

11 

84 

K2 

82 

56 

34 

.9 

169 

58 

00 

7 

8 

224 

9 

323 

2 

20 

. 2 

11 

84 

K3 

82 

49 

26 

.6 

171 

09 

52 

13 

4 

239 

10 

313 

2 

25 


11 

84 

J1 

83 

35 

26. 

.7 

171 

35 

34 

6 

11 

387 

9 

310 

1 

20. 

.4 

11 

84 

J2 

83 

19 

3. 

.8 

173 

04 

26 

3 

15 

333 

11 

314 

2 

20 


11 

84 

J3 

83 

7 

22. 

.4 

174 

55 

29 

18 

16 

385 

6 

329 

1 

27. 

9 

11 

84 

C4 

84 

57 

52. 

.1 

165 

38 

44 

1 

14 

75 

18 

25 

13 

23. 

0 

11 

84 

G3 

83 

25 

42. 

,3 

162 

43 

12 

5 

9 

312 

10 

338 

2 

21 


11 

84 

G4 

83 

16 

40. 

.0 

161 

26 

55 

2 

20 

264 

14 

342 

3 

23. 

2 

11 

84 

LI 

83 

06 

11. 

8 

172 

25 

21 

25 

40 

15 

5 

333 

17 

26. 

2 

11 

84 

LP1 

83 

00 

04. 

3 

172 

59 

57 

9 

20 

16 

7 

126 

24 

28. 

1 

11 

84 

0 

83 

47 

14. 

0 

166 

01 

28 

8 

10 

369 

8 

304 

1 

30. 

9 

11 

84 

CAMP 

83 

37 

14. 

4 

166 

44 

31 

17 

28 

245 

5 

318 

1 

17. 

0 

11 

84 


Downstream B 


MO 

84 

17 

45 

.9 

158 

10 

58 

12 

23 

471 

6 

285 

1 


4 


12 

84 

DNB 

84 

10 

27 

.6 

154 

18 

43 

45 

23 

517 

5 

269 

.5 0. 

.5 

19 

.4 

12 

84 

E19 

84 

9 

3 

. 1 

156 

51 

07 

5 

38 

465 

8 

279 

1 


7 

. 2 

12 

84 

G1 

84 

3 

57 

.6 

152 

08 

59 

5 

4 

534 

12 

269 

1 


12. 

. 3 

12 

84 

G2 

84 

0 

28 

.2 

150 

32 

03 

17 

11 

551 

7 

259. 

9 0. 

.7 

15. 

. 9 

12 

84 

A19 

84 

24 

26. 

. 5 

154 

42 

21 

15 

20 

487 

6 

270. 

6 0. 

. 7 

11. 

2 

12 

84 

B18 

83 

57 

00. 

.1 

153 

33 

32 

12 

23 

426 

6 

279. 

.1 0. 

8 

9. 

0 

12 

84 

B25 

83 

52 

01. 

.8 

153 

10 

38 

4 

23 

13 

9 

268 

38 


9. 

3 

12 

84 

H2 

83 

53 

49. 

.0 

150 

25 

25 

19 

38 

30 

5 

203 

10 


17. 

3 

12 

84 

Downstream 

C 

















H5 

82 

35 

25 


153 

14 

54 


27 

6 

.5 

321 

7 


12 


12 

84 


All azimuths are measured clockwise from true North 
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TABLE 2 VELOCITIES ACROSS ICE STREAM B MEASURED BY TRANSECT 


STATN LAT LONG X 



deg/mi 

n/s 

ec 

deg/min/sec 

km 

B19 

83 

55 

58. 

1 

153 

24 

01 

301. 

68 

B18 

83 

56 

56. 

9 

153 

30 

11 

299. 

79 

B17 

83 

57 

38. 

0 

153 

34 

07 

298. 

54 

B16 

83 

58 

18. 

6 

153 

38 

04 

297 . 

29 

B15 

83 

58 

59. 

9 

153 

42 

06 

296. 

02 

B14 

83 

59 

42. 

8 

153 

46 

20 

294. 

69 

B13 

84 

00 

21. 

9 

153 

50 

13 

293. 

49 

B12 

84 

01 

03. 

9 

153 

54 

24 

292. 

19 

Bll 

84 

01 

45. 

6 

153 

58 

33 

290. 

90 

BIO 

84 

02 

42. 

2 

154 

04 

15 

289. 

15 

B9 

84 

03 

29. 

,5 

154 

05 

25 

288. 

31 

B8 

84 

04 

14. 

9 

154 

06 

32 

287 . 

50 

B7 

84 

04 

59. 

,8 

154 

07 

40 

286. 

70 

B6 

84 

05 

46 . 

.0 

154 

08 

49 

285. 

88 

B5 

84 

06 

32. 

.9 

154 

10 

00 

285. 

05 

B4 

84 

07 

19. 

.8 

154 

11 

11 

284. 

21 

B3 

84 

08 

06. 

.2 

154 

12 

22 

283. 

.39 

B2 

84 

08 

52 

.5 

154 

13 

33 

282. 

.56 

B1 

84 

09 

41 

.2 

154 

14 

48 

281. 

.70 

DNB 








280 

.87 

Al 

84 

11 

10 

. 3 

154 

17 

06 

280 

. 11 

A2 

84 

11 

56 

. 7 

154 

18 

19 

279 

.28 

A3 

84 

12 

41 

.7 

154 

19 

28 

278 

.49 

A4 ; 

84 

13 

25 

.6 

154 

20 

37 

277 

.71 

A5 

84 

14 

12 

.3 

154 

21 

50 

276 

.88 

A6 

84 

14 

45 

.4 

154 

22 

43 

276 

.29 

a7 

84 

15 

31 

.9 

154 

23 

56 

275 

.46 

A8 

84 

16 

02 

.2 

154 

24 

44 

274 

.92 

A9 

84 

16 

49 

.0 

154 

26 

04 

274 

.08 

A10 

84 

17 

19 

. 3 

154 

26 

57 

273 

. 53 

All 

84 

18 

05 

.1 

154 

28 

15 

272 

.70 


Y VELOCITY 



Magnitude 

kz im . 

Bearing 

km 

m/a 


deg 

deg X 

602.45 

386 . 1+ 

0.6 

274.3 

210.9 

601.36 

431.2 

0.5 

274.3 

210.8 

600.57 

445.9 

0.5 

274.3 

210.7 

599.79 

455.4 

0.4 

274.1 

210.5 

598.99 

462.3 

0.4 

274.0 

210.3 

598.17 

467.6 

0.4 

273.9 

210.1 

597.41 

471.8 

0.4 

273.6 

209.8 

596.61 

475.9 

0.3 

273.5 

209.6 

595.80 

480.0 

0.3 

273.3 

209.3 

594.71 

485.6 

0.3 

273.1 

210.0 

593.49 

490.4 

0.2 

272.7 

208.6 

592.32 

494.7 

0.2 

272.5 

208.3 

591.17 

498.8 

0.2 

272.0 

207.9 

589.98 

502.2 

0.2 

271.6 

207 . 5 

588.78 

506.2 

0.2 

271.2 

207.0 

587.57 

509.8 

0.1 

270.7 

206 . 5 

586.38 

512.9 

0.1 

270.3 

206.1 

585.19 

515.0 

0.1 

270.0 

205.8 

583.93 

516.2 

0.1 

269.8 

205.5 

582.74 

517.0 

0.1 

269.0 


581.64 

517.7 

0.1 

269 . 6 

205 . 3 

580.45 

518.0 

0.1 

269.4 

205.1 

579.29 

517.8 

0.1 

269.0 

204.7 

578.16 

518.1 

0.2 

268.7 

204.3 

576.96 

517.4 

0.2 

268.3 

204.0 

576.10 

517.1 

0.2 

268.1 

203.7 

574.91 

516.2 

0.2 

267.9 

203 . 5 

574.18 

515.7 

0.2 

267 . 8 

203.4 

572.93 

516.0 

0.3 

267.7 

203.3 

572.15 

516.8 

0.3 

267 . 6 

203.2 

570.98 

517.7 

0.3 

267.7 

203 . 3 


^Bearing is given as angle clockwise from positive X-axis. 
At station Al , positive X-axis is oriented 64.3 degrees 
clockwise from true North. 
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Figure 3: Velocities across Ice Stream B. 



DISTANCE FROM B25 (km) 



SECTION 2: STRAIN RATE AND ROTATION 


The rate of deformation at the surface of an ice sheet with 
respect to- mutually perpendicular axes (x, y and z) can be 
expressed by the sum of two second-rank tensors, assuming there 
is no rotation in the vertical plane, 



where the strain rate tensor is given by 

Eij = 0.5 ( dUj^ + dUj ) 
dj di 


and the vorticity tensor is 

= 0.5 ( dl^ - dUj ) 

dj di 


and Ui is the velocity in the i direction. The vertical strain 
rate is derived from the incompressibility of ice and the sum of 
the two principal surface strain rates: E zz ~ “ ( E xx + E yy) • 

To obtain strain rates in an area, two methods were employed. 
Either strain rosettes or double lines of stakes were set out 
and surveyed twice; the second survey typically followed the 
initial survey after a year. In the case of strain rosettes, the 
strain rate is determined by assuming a linear velocity field 
between any three stakes in the rosette (usually the three 
outlying stakes) . Velocities are measured relative to the central 



stake, and are given by the horizontal change in position of the 
outer stakes divided by the time interval between surveys. 

The velocity within the rosette is then expressed in the form, 

U x = Alx + A2y + A3 

U y = Blx + B2y + B3 

where x and y are position coordinates and U x and U y are the 
velocity components in the x and y directions. The position and 
velocities are known for the three outlying stakes so the 
constants A1 to B3 can be determined from solving a set of 
simultaneous equations in the form. 



where the letters in parentheses denote the three outlying 
stakes. There is a similar expression for the components in the y 
direction. Assuming infinitesimal strain, the strain rates are 
then given by: 


A1 = dU x /dx = E xx 
B2 = dU y /dy = E yy 

0 . 5 ( A2 + Bl) = 0.5(dU /dy + dU v /dx) = E v „ 

y xy 

where positive values denote extension and negative values denote 
compression. 

For the double line of stakes, the first stage in calculating the 
strain rates is again to determine the velocity field. This is 
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done by solving an overdetermined set of time-dependent 
observation equations to derive, with, error estimates, the 
position and velocity for each stake site, using the method 
described in Wager et al (1980) . The position and velocity of the 
stakes are then used in a series of interconnecting triangles 
down the scheme using a different linear velocity field for each 
triangle. 

Measurement Errors 

A combination of field measurements of distances, angles and 
azimuths was made to determine strain rates. During the first 
field season, 1983-84, distances were measured by CA 1000 
Telurometers and HP Electronic Distance Meters (EDM's), which are 
accurate to within 3 cm per 1500 m and 1 cm per 1500 m, 
respectively. Geodimeter 112 EDM/s were used during the next two 
field seasons; accurate to _+ 1 cm per 1500 m. The greater 

accuracy of the EDM's made their horizontal positions above the 
stake significant. The theodolite mounting of the EDM was 
optically plumbed onto the stake, so we consider a reasonable 
error to be 1 cm. However, where strain nets were measured using 
only distances the misclosure was more typically 2 cm. 

The mean misclosure between rounds of the angle measurements was 
4 seconds of arc. However, one of the instruments in the 1984-85 
season consistently misclosed to 20 or 30 seconds of arc. In this 
case, measurements were repeated up to five times to obtain a 
consistent reading. In calculating the errors of the strain rate, 
a minimum error of 5 seconds of arc was used; if the misclosure 
was greater than 5 seconds, the size of the misclosure was taken 
as the value of the error. 

Azimuths were measured with the least accuracy. The greatest 
contributions to the azimuth error were in the determinations of 
time and longitude. A 1-second error in time contributes to about 
a 15-second-of-arc error in azimuth; a l-second error in 
longitude gives a 1-second-of-arc error in azimuth. Time was 
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taken from the satellite tracking unit, or from radio time 
signals. Sometimes wristwatches were used which were later 
calibrated using a satellite tracking unit. We have not 
calibrated the time given by these units, but we note that when 
they have been run concurrently they can differ in time by up to 
1 second. Therefore, we have used an error of 1 second for time. 

The position in longitude is known to better than 1 minute when a 
geoceiver position is used, but when single SATNAV fixes are used 
the longitude error could be up to 10 minutes of arc. SATNAV 
satellite receivers (Racal-Decca 412) are a single -channel system 
accurate to ^200 m after about 5 passes, and were used in the 
place of the MX1502 geoceivers during the 1983-84 field season. 
Only a few station positions have been estimated by one SATNAV 
position; the remainder are determined using either three SATNAV 
fixes or up to 50 geoceiver fixes. We have used 45 seconds of arc 
as the error in azimuth. This may be high for azimuths measured 
under optimum conditions: where time is taken directly from a 
geoceiver using a stopwatch to measure lapse time, and longitude 
is taken from a geoceiver position. Figure 4 and Table 3 present 
the strain rate data. 
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Figure 4: Regional map showing selected strain rates. 


L5 



*3 

ro 

ro 

0 

N 

r~> 

rr 

II 

- 

P 


ro 

rr 

1 

ro 

H> 

,-~N 

II 

O 

X 


P 

r— * 

ro 


X 

rf 


X 

M* 

to 

CO 

a 



O 

0 

- 

M- 

h-* 


X) 

O 

< 

P 

O 

<D 

I— 1 

7T 

»T 


£ 

rr 

P 

H- 

1 — * ■ 

X 

Co 

0 

a> 

0) 

P 

CO 

- 

r— 1 




0 

P 

10 

My 

D 

rr 


a 

r< 

zr 


p 

0 

cr 

r-*- 

rt 

fD 

P 

M’ 

p 


N 

n 

rr 

O 

H* 

p 

P 

D 

rr 

rr 

CTQ 

fD 

P 



r — * 

M> 



(0 


rr 



(V 

O 


p 

r— • 


CO 

O 


r j. 

n 


O 

?r 


P 

*: 



m* 


P 

co 


P 

CD 


0- 

rh 


n 

r* 


0 

O 


g 

3 


ro 



>r 

rr 


0) 

rj 


w 

C 


trt 

fD 


f“*' 



0 

Z 


p 

0 



r{ 



rr 


n 

:r 


0 


4 >CO^^Cn*>coMh^OZZf 0 ^l^^^rOH^t^£i~* 

W CO ^OOL/'OOOvO^DH 


10 

X) 

fD 

n 

rr 


m 
H 
50 

Ocnp'OOOOWOOrOHH^WhOHMOMfONJO^ 

^r^roc/icncocncoco^OvOr^co\^cocriro^cj\^cnv£>no 

^H^hOrONivDWONOOUlOsJWUiHWOOWO^O 

OOOOOOOOOOOOOOOOOOOOOOO 

ooooooooooooooooomooooo 

cn^c^Nocofo^^h^rocOc/iNo^^rocooforocorofo 


n 

00r^7s^^c H c^xiQ0n00 ,p ii^ T jp]r t it 1 it I irntT]piPiDDnnnnhflmi. f i>^ 
^N^oroi-rot-^co^ror-roM ^£S£SS£S 22222 S££££ 

cnw 

^ ;N 0 ^ WH 000 ^^ 0 mijwui ^^ u 1 ^ 00 v 0 vD ^^hwOhhwhwu;|> 

^^IOC/i^h^4>^OUiONvOC04>C*OCo4>LnvO<»CC^^cO^^O?Co2cri)2vog 

00 000000 000 000000 °000000000000000000 
^^° 2 °!^ 00000000000<0<: ^ >< ^ :>00c=::i OOOOOOOOOOOOO 

4>^roh^^rorouirororsjvocoi-*rorororouJLnoo(-^rounfsjp-N3u>(NororoSSStsj^ 




C/> 


MlMjiC^ro^Ol — ‘ ( — , N)OvDOOh- lN vJtJvjD| — > Q CO £> I> 
NJCJ^^>OOOJ^jOOhO-^jO OACC^COrOfOOC^CncO'OrO 

00000000000000000000000 

ooooooooooHooooooooooob 

O^LM^rocnrocriroi — > ro co^roro.f>cocorococo £> ^ 


I II 4 14 11 

OOwOOOOOOOOOOroMOOro|>ON)NO 

UJWxiMKJH^^i^^^oWHNJvOfOUi^CXJOO^ 
4>C/iCOCTNcO.(>vI}vOr— 00 ONONirOvOvCOUlOOHOOOO 

OOOOOOOOOOOOOOOOOOOOOOO 

oooooooooo^oooooooooooo 

\j^oomuimcm*jmuj a>^horo^.p-cocococo.p'oo 


O’sJOOl-'HOHMOOWHl-'WCJhOC^OCDUlHl-'M 

h-vONJ^^^OvlCTU^MOQLnM^O^vDC^OUi^ 

OOOOOOOOOOOOOOOOOOOO 

rororji^roNJi^roroioioroiorororororoioro 


C0^^^CXDC04>.p'CO^c^a>00'O00vO.f>c0C/ic0C^a>tCrs 
CnC/lC0roa>C0Or^^h^^C/lC/l.J>00L0c0'0C04>»— ‘lOfO 

r— * 00 co ro ^>^00 Oi — , u)^h- i Ni|>wo^)oo^ja> 

VO 

O' ‘i ‘OOOi — ‘OOOOOOOOOOOOO 

oc ' * ro \c co co r — » (— * co co .p" 0 ro co co o> 

co 


0 OWW 0 Hww ^WUi^^OOHNO>^MWOOONJOMOro 
S^^!^ WC0 ^ O ^ aOOOCD ^^^WOOCX)H(^a\v0v0OCD(»0>v0U>cnvr) 

UiO00NiWWWcn^(»UJWu^O^W^vOCncn00\r)00^t0C» W|sjxja>s^co^W 

I + 

° ° ° ° ° ° 0 0 0 0 0 0 0 0 0 0 ° ° O O O O O O O O O O O o O O O O O O 

O00000000000000000000000000000r^r^r~ir-!>r^r^ 

p-p-roco^ioroc/i[ororoar' l roh-foroiorococna>cjroa'Nro.(>[oro[ororo(orororoh^ 


OOOWh ‘OOOOr-^OOOOOOOOOO 


OOror^i-^OOOOOOOOOco^ 


I + 

00 000000 000000 000 ° 00 0000 000000000000 

OOOOOOOOOOOh-*OOOOOOOOh- ‘OOOOOOOOor^r^r^r^r^r' 

Chcncoco^roco^coroco NoroNocococop>2S2SgS22£2S22222S 


HH °!^,^ HW f W0 [° HMr0W ^ HH ° M O0U^WHHCni!nOHOOHh 1 O 
C0IJlh-W^vJO^r0Wv£)00WOWUlC0^UlWON)00OOCD^HM^M0>vJCn^^ 

| + 

000000000000000000 ° 00 000000000000000 

roro^roioro,ororoioio^^ioroioro^NjiohoioN3(o(ofoioiorororofororororo 


^ 4>4>aN^p'vO^^C/iCncncnOGCn4>cn4>P>c/iP'' 

wwvONjwwcocnCDONOro^MMUiw^OO^OO^HOw^HO^^fotoi' 


CO Ul 


WNlO 


cr* co 


MNjvirosi^oowvo 

O ro ro 00 h-* 


O'* cr\ -o o\ o P> 

vD h-» 


ctn 4> O ^ ro 4> cO 

00 

oowoooHHOoqrooooooooouiroooooooMOoooocjo 

r^^ro Uir-H coooOMOOhO oi-oJ-fOK- k> Ij 2 i -- 2 O 

cn-ocn -vjvO ro p- C- 


C : 

►-i 


O 


X 
h— * 

o 


H 

> 

X DO 

5 K 

. ra 

Co fO U) 
P * •.* 

i 

I-* W 

H 


5 

M 

Z 

2 


H 

* x) m 

WN CO 


X ?0 

r- o 
O H 
• > 
co H 

ri 
P 

a 


p 


o 

»-n 


o 

z 

+ 


CO 


OJ CTj CXy a; Cc CO CO (X a Co 03 OO CO CO 00 Co Cij 00 CD CD CO CO 
U1 ^ ^ *' *' t' ^ ^ co ^ i > 4 > coco 4 >^ co ? S 2 

cocococococccqcoccoooocooooooogooooooooo 

c/i cm cm c/i c/i lyi c/i c/i ui cn c/i 4 > u » cn |> cn i > p -. 


“ 2 ““ l - SSS ? S 22 SSP ^^, OOOOCOOO ( 5?? c,:>coa)00 00 00coooc,o « ,co «“<-'c 


S ? S g S X g “ 5 £ £ S*SS "? 55 *“*®"“"® S 


O 

> 

H 

r>: 





SECTION 3: SURFACE ELEVATIONS FROM GEOCEIVER AND OPTICAL LEVELLING 


The Doppler satellite tracking system also provided measurements 
of surface elevation with respect to the reference ellipsoid WGS 
72. This ellipsoid has a semi-major axis of 6378135 m and a 
flattening of 1/298.26 . Table 4 gives the measured elevation 
with each position determination. At 24 stations the elevation 
was measured twice; in 21 cases, the agreement was 8 m or less. 
Again, the standard deviation of elevation provided by the MX1502 
decreased as the number of passes increased in the same manner as 
for the horizontal position. Seventy-five percent of the 
differences were within one standard deviation of the least 
accurate of the two elevations, and no differences were more than 
two standard deviations of the least accurate elevation. 

Optical levelling was also carried out to gain detailed profiles 
of the ice surface. The distance between the level rod and 
optical level varied from 150 to 250 m. In early surveys, 
distances were estimated by Skidoo odometer; some were remeasured 
using EDM's as part of the strain-rate survey. During later 
surveys, upper and lower cross hairs were read which gave the 
distance between level rod positions and also served as a check 
on errors. Level lines were not closed. On three occasions, 
optical levelling was carried out between stations which had 
elevations that were also determined by geoceiver. The difference 
in elevation determined by each method typically misclosed by 2 
to 3 m. This misclosure is fully accounted for by the error in 
the geoceiver-determined height and is consistent with the given 
errors. Table 5 gives geographic position of each profile; Table 
6 gives the data. Figure 5 shows the geographic locations of each 
levelling line and Figure 6 presents the elevation profiles for 
these lines. 


17 



c; 

LO ui os vO l/l 
O 
O 
P 
w 

p 


vOO I—* 


a 

XDnncsDo>>?*DO 

UiHU>HWN)UN)|“‘2;i 
VO op 

w 
rr 
M 
CD 
CD 
3 


a 

^wrOHWhOMroH '(-“•toi-'S; t— 1 

LncovD CO v£) p 


(n 

rr 

0) 

CD 

3 


n 


MOOOOUK-'vOOO'vIDO 

Ln0J^>Oui^>h-*004> 


UHudNJMN)Ov^OOvO\£) 

OOUUMOj>C>OVHvD 


^^ai^Ln-vjH^oooouJUJvooooo 

O'-JfOvO0 , >O00CT''^M^U)OUlh J 


o 

xonooorroon^^^^^^nnn 

ro ro oj ro ro rohjH-‘£'OJ4>oj{OM ‘ wmh 1 -C>locu 

. ^ m 

> 


«yi 

rr 

CD 

a 

h*« 

O 

a 


to -F- vo ro io - o uj ro cj >— * ro crv ut 

O^OOQOUIvO^J W^HOOOOC^POO 


io 

t_n oo ui 4> -C> £> p'vou , ><oL->h- i rooco i— ‘roCT'P" 

OvO(^UiCOPLnP^WWK)^(-ivjpc»OavUi 



i—* 4> h-* 

v£> PfOUI'vJLnUll/lK) 


i — * ro h-» m * h- 1 

CX3^vOLnrOtOh^OOOOvU>'OUiLncrv 


3 

h- 1 


z: 

o' 


a> 


co 


ro roro 
Ou^wro 

4> h- i, 'OHH roro 

OouiOoo^ouiuiCTv.p^on 


roro to Loiororoh- " rotoro 

ooooooP CJOWWOl-'PWOOU 


MhP uu— * rop"fo f — * i — ‘ r— » i — 1 oj 

LnrOvlvJCDOvOOOOvDPOVUiHPOOOOvJO 


$ 


o 

f-h 

T3 

CD 

w 

w 

fl> 

w 


a 


sj vJWMHvOl — O nI 


to h- 

oovyiX>h-*cx34>P v oDLncrP'OCrvrooj 





(D 

P 

O 

<D 


ui o av O'. 


o (o oi oo a> 


OOO-vJLnPOOO^HVOO^OCOCO 


m 

m 

o 


n 



a 

m* 

r-n 

M-* Ml 

a> m 

M 

rt 

o 

n 


TABLE 4: ELEVATIONS FROM GEOCEIVER MEASUREMENTS 
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TABLE 5: POSITIONS OF OPTICALLY LEVELLED TRAVERSES 
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RADIO-ECHO SOUNDING OF ICE THICKNESS 


SECTION 4: 


Spot measurements of ice thickness were made at several survey 
stations (Table 7) and along three short profiles near the 
grounding line of Ice Stream B (Tables 8, 9 and 10 and Figure 7). 

The radio-echo sounder used was on loan from the University of 
Wisconsin, Geophysics and Polar Science Center. A description of 
this radio-echo system is given by Shabtaie and Bentley (1987) 

To reduce the two-way travel time to ice thickness, a velocity of 
169 xlO' 6 m/s was used for the radar energy (Shabtaie, personal 
communication) . 
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TABLE 9: THICKNESS PROFILE FROM A2 (CIR) TOWARD N1 












SECTION 5 : ACCUMULATION RATE 


Two methods were used to nisasure the accumulation rate. The 
preferred method involved coring to a depth of about 10 m and 
sampling the core for beta particle activity. There are peaks in 
this measurement in the 1956 and 1965 horizons associated with 
global nuclear bomb tests. Cores have been taken at four sites, 
(DNB, CIR, E19 and DNC) but the analysis has not been completed. 

The second method is to repeat measurements of the length of 
exposed survey stakes. The normal time interval between 
measurements of stake heights was only 1 year, so the measurement 
of accumulation rate may not be representative of the longer term 
average. To convert these stake exposure measurements to mass 
equivalents, we suggest a value of 388 kg/m -3 as the average 
density of the surface firn. 


Table 11 includes these measurements indicating location and the 
time interval. The standard deviation of the accumulation rate is 
derived from the eight stake exposures (four survey poles and 
four adjacent bamboo poles) measured at each strain rosette 
covering an area of typically 7 km 2 . The stake lines at 
Downstream B usually consist of 10 stakes along an 15 -km line. 


TABLE 11: ACCUMULATION RATES FROM STAKE EXPOSURES AT STRAIN ROSETTES 
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SECTION 6: TEN-METER TEMPERATURE AND FIRN DENSITY PROFILES 


At the locations where 10-meter-long cores were removed, 
temperature measurements were made in the resulting holes. In 
1983, a single thermistor was used at a single depth, while in 
1984 and 1985 a thermistor chain was hung in the hole to measure 
the vertical temperature gradient. The single thermistors were 
glass bead Fenwall GB32M2 (2000- ohm resistance at 25 degrees C) . 
The thermistor chain alternated these with Fenwall GB41P2 
thermistors (10,000-ohm resistance at 25 degrees C) every meter 
for 5 m. Each thermistor was calibrated at zero and three sub- 
zero temperatures. A precision Wheatston bridge (Leeds and 
Northrup Model 4289-3) was used. The hole was covered to reduce 
air circulation. Repeated measurements were made as soon as the 
core was removed to allow the cooling curve to be determined, 
giving a more accurate estimate of the equilibrium temperature. 
Measurements were continued until the cooling rate was less than 
0.01 degrees C per hour, which usually occurred after 4 to 6 
hours. When possible, a final measurement was made after many 
hours to confirm the calculated equilibrium temperature. Self- 
heating of the thermistors caused by the 1-mA current output by 
the measurement bridge limited the accuracy of these temperature 
measurements to_^_0.05 degrees C. Table 12 presents these data. 

Firn density profiles were calculated using volume and weight 
data collected at sites where 10-meter-long cores were extracted 
(Figure 8 and Table 13) . In the process of preparing the 10-meter 
core for shipment, each core was cut into 10— cm sections and 
weighed to the nearest tenth of a gram using a sheltered triple- 
beam balance. During the 1983 and 1984 field seasons the 
weighing was done after each 10-cm section had been placed in a 
small plastic bag, while in 1985 the core sections were weighed 
beforehand. The weight of the bags is considered negligible and 
within the error of these measurements; hence no attempt is made 
here to adjust the 1983 and 1984 data. In 1983 and 1984 the 
diameter of the core was taken to be that of the inside diameter 
of the coring assembly (7.3 cm) , while in the 1985 season, each 
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core section was measured individually using a caliper. The 1985 
data showed that the core's diameter varied from the inside 
diameter of the coring assembly by less than 2mm. Density values 
obtained from the 10-cm segments of core are assigned to the 
vertical position of the center point of each segment. 
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TABLE 12: TEMPERATURE MEASUREMENTS 
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Firn Density Profile At CIR Camp (Dec. 1984) 



.32 .34 .36 .38 .4 .42 .44 .46 .48 .5 


.3 


.52 .54 .56 


Density (gm cm-3) 


Figure 8 (a-c): Firn density profiles. 
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Firn Density Profile At DNB Camp (Dec. 1984) 
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Firn Density Profile At E19 (DNB) Camp (Dec. 1985) 
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SECTION 7 : SPECIAL SITES 
RESULTS FROM DOUBLE STAKE LINES AT LI AND DNB 


Special sites refer to two areas where double stake lines were 
established to determine the detailed strain field in an extended 
region. The relative velocities and strain rate were determined 
as previously described in Section 2 
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TABLE 14: VELOCITY PROFILE FROM DOUBLE LINE OF STAKES WHICH CROSSES 
THE EDGE OF CRARY ICE RISE NEAR STATION LI (see Fig 11(a)) 


STATION 

COORDINATES* 

VELOCITY 

NAME 

X 

Y 

X-comp . 

Y- comp 


(m) 

(m) 

(m/a) 

(m/a) 

W 

50000.00 

50000.00 

0.00 

0.00 

X 

49590.38 

50187.44 

0.51 

1.06 

A3 

50468.18 

50180.27 

0.52 

-1.13 

Y 

49741.70 

50964.61 

2.75 

1.02 

A2 

50505.86 

50981.21 

2.73 

-0.88 

Z 

49827 . 13 

51614.92 

4.78 

1.22 

AL 

50536.26 

51583.82 

4.67 

-0.63 

AO 

49860.28 

52300.72 

9.66 

1.46 

0 

50572.51 

52265.96 

8.56 

-0.41 

A 

49879.28 

52737.89 

22.90 

2.33 

L 

50782.65 

52820.14 

23.85 

-0.02 

B 

49886.65 

53322.97 

42.55 

3.56 

2 

50752.39 

53544.42 

48.81 

1.67 

C 

49914.35 

54282.94 

74.44 

5.04 

3 

50747.48 

54509 . 57 

81.04 

2.21 

D 

49954.55 

55430.03 

109.08 

5.00 

4 

50852.28 

55476 . 13 

110.58 

2.31 

E 

49991.85 

56612 . 34 

141.96 

4.56 

3 

50989.94 

56690.43 

144.82 

1.13 

F 

50027.91 

57774.16 

173.96 

3.52 

6 

51043.10 

57828.09 

176.34 

0 . 46 

G 

50062.47 

58725.45 

197.69 

2 . 64 

7 

51085.83 

58867.63 

202.36 

-0.56 

H 

50093.00 

59677.40 

219.57 

1.96 

8 

51135 . 38 

59928.41 

226.66 

-1.83 

I 

50125.66 

60685 . 17 

241.91 

1.02 

9 

51175.01 

60694.78 

243.77 

-2.77 

J 

50190.23 

61808.99 

265 . 53 

-0.33 

10 

51231.68 

61743.33 

265.92 

-4.23 


*Aziriiuth of X-axis is 316.9 degrees. Geographic coor- 
dinates of station W are: S83°06'll", W172° 25'23". 



TABLE 15: STRAIN RATE PROFILE FROM DOUBLE LINE OF STAKES WHICH CROSSES 
THE EDGE OF CRARY ICE RISE AT STATION LI (see Fig 9) 


STATIONS 

COORDINATES* 

PI 

f 

po + 


P2 

f 

ANGLE OF 

PI 




X(m) 

Y(m) 

(io- 

V 1 ) 

(10 a 

r 1 ) 

(10* 4 

a' 1 ) 

WRT X - ci 

. X j 

.. .s 

W 

X 

A3 

50019 

50122 

0*22: 

0.2 

-0.0 +0.2 

-0.2- 

0. 3 

59 . 3 + 

0 

1 

J 

X 

A3 

Y 

49933 

50444 

0.52 

0.06 

-0.04 

0.09 

-0.5 

0.1 

67.74 

0 

08 

A3 

Y 

A 2 

50238 

50708 

0,45 

0.06 

-0.12 

0.09 

-0.3 

0.1 

76. 70 

0 

08 

Y 

A 2 

Z 

50025 

51187 

0.75 

0.07 

-0.2 

0.1 

-0.5 

0.1 

69.65 

0 

05 

A 2 

Z 

Al 

50290 

51393 

0.69 

0.08 

-0.1 

0.1 

-0.5 

0.1 

65 . 18 

0 

07 

Z 

A1 

AO 

50074 

51833 

2.59 

0.08 

-1.95 

0.08 

-0.6 

0.1 

46 . 99 

0 

01 

A1 

AO 

0 

50323 

52050 

1.39 

0.09 

-2.21 

0.08 

0.8 

0.1 

59.48 

0 

01 

AO 

0 

A 

50104 

52434 

14.9 

0.1 

-12.9 

0.1 

-2.0 

0.1 

4 7 . 24 

0 

00 

0 

A 

1 

50411 

52608 

12.92 

0.08 

-12.69 

0.08 

-0.2 

0.1 

48 .65 

0 

.00 

A 

1 

B 

50182 

52960 

15.61 

0.08 

-15.49 

0.08 

-0.1 

0.1 

48 . 83 

0 

.00 

1 

B 

2 

50473 

53229 

16.26 

0.07 

-15.60 

0.07 

-0.65 

0.09 

48.41 

0 

00 

B 

2 

C 

50184 

53716 

15.57 

0.06 

-15.23 

0.06 

-0.34 

0.08 

47 . 70 

0 

00 

2 

C 

3 

50471 

54112 

14.63 

0.06 

-15.26 

0.06 

0.63 

0.08 

46 . 64 

0 

00 

C 

3 

D 

50205 

54704 

13.30 

0.06 

-13.50 

0.05 

0.21 

0.08 

45.42 

0 

00 

3 

D 

4 

50518 

55138 

14.04 

0.06 

-13.50 

0.06 

-0.54 

0.08 

4 5.35 

0 

00 

D 

4 

E 

50266 

55839 

12.39 

0.05 

-12.43 

0.05 

0.04 

0.08 

44 . 40 

0 

00 

4 

E 

3 

50611 

56259 

12.42 

0.05 

-12.34 

0.05 

-0.07 

0.07 

4 3 56 

0 

00 

E 

5 

F 

50336 

57025 

12.05 

0.05 

-12.14 

0.05 

0.86 

0.07 

43.23 

0 . 

00 

3 

F 

6 

50687 

57430 

12.57 

0.05 

-12.14 

0.05 

-0.43 

0.07 

43.45 

0 . 

00 

V 

l 

6 

G 

50377 

58109 

11.12 

0.05 

-10.90 

0.06 

-0.21 

0 . 08 

42 . 60 

o 

. Oi) 

6 

G 

7 

50730 

58474 

11.14 

0.05 

-10.91 

0.05 

-0.22 

0.07 

4 2 .46 

0 

00 

G 

7 

H 

50413 

59090 

10.37 

0.06 

-9.62 

0.06 

-0. 76 

0.08 

42 . 15 

0. 

00 

1 

H 

8 

50771 

59491 

9 . 92 

0.05 

-9.67 

0.05 

-0.25 

0.07 

41 . 59 

0 , 

00 

H 

8 

I . 

50451 

60097 

9.72 

0.05 

-9.09 

0.05 

- 0 . 64 

0.08 

41 . 50 

0 

00 

8 

I 

9 

50812 

60436 

9.68 

0.06 

-9.15 

0.07 

-0.53 

0.09 

41.03 

0 

00 

I 

9 

J 

50497 

61063 

9.05 

0.06 

-8.47 

0.07 

-0.58 

0.07 

40 . 7 S 

0 

00 

9 

J 

10 

50865 

61416 

8.98 

0.07 

-8.47 

0.06 

-0.50 

0.08 

40.23 

0 

00 


"Azimuth of X-axis is 316.9 degrees 

+ P1 and P2 are principal axes of horizontal tension and 

compression, respectively. P3= -(P1+P2) is vertical strain rate 
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LI Scheme 



Figure 9: Strain rates measured along a double stake line at Li, across the margin of Crary Ice Rise 


TABLE L6 : VELOCITY PROFILE FROM DOUBLE LINE OF STAKES PARALLEL 

TO FLOW IN THE MOUTH OF ICE STREAM B. 


STATION COORDINATES* VELOCITY 


NAME 

X 

(m) 

Y 

(m) 

X-ccmp . 
(m/a) 

Y-comp 

(m/a) 

F10 

48722.25 

65058.83 

- 47.29 

495.84 

E10 

47123.69 

64442 . 66 

-53.71 

498.97 

F9 

48667 . 88 

63612 . 12 

-51.58 

497.55 

E9 

47083.31 

63002.17 

-57.15 

500.27 

F8 

48614 . 64 

62199.78 

-55.19 

498.63 

E8 

47043 . 61 

61580.14 

-60.07 

501.56 

F7 

48559 . 75 

60746.05 

-58.74 

499.84 

E7 

47004.03 

60149.06 

-63.26 

502.21 

F6 

48504.94 

59312.11 

-61.99 

501.06 

E6 

46965.05 

58715.93 

-65 . 90 

502.75 

F5 

48450 . 11 

57894.97 

-64.67 

502.02 

E5 

46925.04 

57248.23 

-68.72 

503.54 

F4 

48394.46 

56468.13 

-67.50 

502.78 

E4 

46885 . 73 

55817.89 

-70.98 

504.50 

F3 

48338 . 80 

55037.41 

-70.06 

504.54 

E3 

46845 . 75 

54354.00 

-74.02 

505.78 

F2 

48283.48 

53615.69 

-73.13 

506.66 

E2 

46804.46 

52842.80 

-78.67 

507.22 

FI 

48238 . 29 

52169.61 

-77.17 

508.23 

El 

46765 . 76 

51429 . 91 

-83.15 

508.76 

A1 

48660.00 

50000.00 

-83.96 

510.92 

A2 

47208.94 

50000.00 

-86.21 

510.74 

D1 

48138.31 

49481.45 

-86.68 

511.58 

Cl 

46786 . 96 

48832.68 

-91.37 

511.56 

D2 

48091 . 78 

48221.00 

-91.18 

512.99 

C2 

46861 . 04 

47383.77 

-96.03 

513.81 

D3 

48040.47 

46856.76 

-95.49 

514.64 

C3 

46936.43 

45934 . 86 

-100.61 

515.35 

D4 

47991.46 

45543.27 

-99 . 90 

516.02 

C4 

47013 . 62 

44481.85 

- 104.74 

516.90 

D5 

47942 .21 

44218.70 

-104.21 

517.35 

C5 

47090 . 69 

43038 . 58 

-108.88 

518.34 

D6 

47892 . 79 

42887 . 81 

-108.17 

518.78 

C6 

47169 . 13 

41569.26 

-113.30 

519.28 

D7 

47843 . 64 

41564.22 

-112.25 

519.62 

C7 

47243.96 

40177 .83 

-118.18 

520.49 

D8 

47778 . 92 

39804.28 

- 118 . 64 

521.08 

C8 

47318.40 

38791.39 

-121.23 

521.17 

D9 

47730.55 

38479 . 92 

-121. 70 

521 . 77 

C9 

47396.86 

37328 . 70 

-125.18 

522.68 

DIO 

47683.85 

37153.12 

-125.62 

523.23 

CIO 

47471.54 

35941 . 58 

-128.17 

523.51 


* Azimuth of X-axis is 289.7 degrees. Geographic coor- 
dinates of station A1 are: S84 ° 10' 52", W154® 10'08 M . 
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TABLE 17: 


STRAIN RATE PROFILE FROM DOUBLE LINE OF STAKES WHICH 
PARALLELS ICE FLOW IN THE MOUTH OF ICE STREAM B (see Fig. 10) 


STATIONS COORDINATES* 
X(m) Y(m) 


F10 

E10 

F9 

48171 

64371 

E10 

F9 

E9 

47625 

63686 

F9 

E9 

F8 

48122 

62938 

E9 

F8 

E8 

47580 

62261 

F8 

E8 

F7 

48073 

61509 

E8 

F7 

E7 

47536 

60825 

F7 

E7 

F6 

48023 

60069 

E7 

F6 

E6 

47491 

59392 

F6 

E6 

F5 

47973 

58641 

E6 

F5 

E5 

47447 

57953 

F5 

E3 

F4 

47923 

57204 

E5 

F4 

E4 

47402 

56511 

F4 

E4 

F3 

47873 

55774 

E4 

F3 

E3 

47357 

55070 

F3 

E3 

F2 

47823 

54336 

E3 

F2 

E2 

47311 

53604 

F2 

E2 

FI 

47775 

52876 

E2 

FI 

El 

47270 

52147 

FI 

El 

Al 

47888 

51200 

El 

Al 

A2 

47545 

50477 

Al 

A2 

D1 

48002 

49827 

A2 

D1 

Cl 

47378 

49438 

D1 

Cl 

D2 

47672 

48845 

Cl 

D2 

C2 

47247 

48146 

D2 

C2 

D3 

47664 

47487 

C2 

D3 

C3 
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C3 
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C3 
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D8 

C8 

47447 

39591 

D8 

C8 
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♦Azimuth of X-axis is 289.7 degrees. Geographic coordinates of 
station A1 are: S84°10'52", W154°10’08". 
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Figure 10: Longitudinal profile of strain rates measured along a double stake line near DNB. 
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SECTION 8: AERIAL PHOTOGRAPHY 


In January 1985 the U. S. Geological Survey took aerial 
photographs of the Crary Ice Rise. The photographs covered an 
area of approximately 175 km by 70 km and were taken at an 
altitude of 25,000 feet. In all, 570 photographs were obtained. 
Using every other photograph, a mosaic was created by matching 
features (stratsugi and crevasses) on adjacent photographs. 
Overlays were then prepared of crevasses, crests of undulations 
and open rifts in the ice shelf downstream of the ice rise. 
Identification of ground control stations provided an estimation 
of the scale to be 1:47,400. The mosaic and overlays were then 
photographically reduced to scales of 1:250,000 and 1:500,000. 
Figures 11(a) and 11(b) show the crevasse and undulation overlays 
with a geographic grid. 

There is a gap in the photographic coverage along a thin strip to 
the grid northeast of the ice rise. Bad weather forced the data 
to be collected in. two missions separated by 3 days. While an 
extra photographic strip was taken to provide ample overlay, poor 
navigation caused errors in position large enough to create the 
gap. Fortunately, there were two control points on the smaller 
mosaic which allowed a separate calculation of scale, and allowed 
this piece to be aligned relative to the larger section. The gap 
is relatively small and its effect on the mapped features is 
small. The most apparent difference in the two data sets are the 
crevasses downstream of the ice rise. In the smaller mosaic, more 
crevasses are apparent compared to the larger mosaic. This is 
most likely a result of differing solar azimuth and elevation at 
the time of each photographic mission. 
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SECTION 9 : MULTI-LEG ROSETTES 


The multi-leg strain rosette technique outlined in MacAyeal 
(1985) was used at three stations near the grounding line of the 
ice stream, Nl, N3 and N4 . A fourth multi— leg rosette was 
planted at N2 , but it was not resurveyed by the end of the 1985 - 
1986 field season. This technique was chosen because it provides 
a means of reducing measurement uncertainty under circumstances 
where the time period between initial station deployment and 
resurvey is short. Our goal was to recover the strain rate data 
after approximately 15 days. A redundancy of 40 measurement legs 
was used by planting 10 outlying stakes and surveying them from 
four independent central stakes. The details of this innovative 
method are discussed in Appendix 3. The data are presented in 
Table 18. Figure 12 shows the similarity of velocities of the 
ten outlying stakes relative to the velocity of four center 
stakes. The computed strain rate is indicated in the center of 
this figure. Figure 13 indicates the scatter of calculated 
strain rates. 
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TABLE 18: MULTI-LEG ROSETTE RESULTS 
L8(a): Comparison between weighted, unweighted and linear method 
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18(b). Comparison of 3-Leg, 10-leg and 40-leg rosette data from N3 + 
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Figure 12-13: Stake movements and strain rates at station N3 from multi-leg rosette results. 
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CAMP 

83 

37 

15 

166 

44 

30 

162 

.484 

<689 

. 595 

6 

12 

33 

1 

27 

47 

g2 


83 

37 

08 

166 

45 

20 

162 

. 366 

-689 

. 845 

6 

12 

41 

1 

27 

43 

g3 

0 

83 

47 

14 

166 

01 

28 

166 

.631 

-669 

.539 

6 

01 

44 

1 

30 

02 

g2 


83 

47 

07 

166 

02 

43 

166 

.440 

-669 

.810 

6 

01 

53 

1 

29 

55 

g3 

A1 

83 

32 

41 

165 

07 

48 

174 

.469 

-657 

.088 

5 

55 

01 

1 

34 

16 

s 1 

A2 

83 

57 

29 

164 

16 

24 

181 

.861 

-645 

.838 

5 

48 

57 

1 

38 

16 

g3 

B1 

83 

35 

13 

165 

58 

09 , 

163 

.688 

-655 

.015 

5 

53 

54 

1 

28 

26 

si 


83 

55 

16 

166 

00 

02 

163 

. 307 

-655 

.015 

5 

53 

54 

1 

28 

14 

s2 

B2 

84 

02 

58 

165 

56 

21 . 

160 

.540 

-640 

.990 

5 

46 

20 

1 

26 

45 

si 


84 

03 

02 

165 

57 

31 

160 

.292 

-640 

.924 

5 

46 

18 

1 

26 

36 

s2 

Cl 

33 

53 

00 

166 

46 

17 

155 

.442 

-661 

. 244 

5 

57 

16 

1 

23 

59 

s2 

C2 

83 

57 

32 

167 

41 

16 

143 

.054 

-655 

.433 

5 

54 

08 

1 

17 

18 

g3 

C3 

84 

26 

08 

166 

57 

08 

139 

! 486 

-601 

.891 

5 

25 

15 

1 

15 

22 ‘ 

g3 

C4 

84 

57 

52 

165 

38 

44 

138 

. 596 

-541 

. 585 

4 

52 

42 

1 

14 

54 

g2 


84 

57 

50 

165 

38 

32 

138 

. 643 

-541 

.637 

4 

52 

44 

1 

14 

56 

g3 

D1 

83 

47 

37 

167 

12 

13 

152 

.660 

-672 

. 133 

6 

03 

08 

1 

22 

c 29 

s 1 


83 

47 

40 

167 

13 

29 

152 

. 392 

-672 

.099 

6 

03 

07 

1 

22 

20 

s2 

D2 

83 

47 

44 

168 

22 

02 

138 

.936 

-674 

. 883 

6 

04 

37 

1 

15 

04 

s 1 

El 

83 

42 

08 

166 

52 

00 

158 

. 921 

-681 

. 128 

6 

07 

59 

1 

25 

52 

r _. i 


83 

42 

25 

166 

52 

29 

158 

. 706 

-680 

. 638 

6 

07 

4 3 

1 

25 

44 

s 2 

E2 

83 

36 

38 

167 

41 

37 

151 

. 249 

-693 

. 317 

6 

14 

34 

1 

21 

43 

si 

El' 

83 

35 

14 

168 

07 

42 

146 

. 518 

-696 

. 986 

6 

16 

32 
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19 

09 

e 1 

E2 ' 

83 

33 

32 

167 

52 

18 

150 

. 301 

-699. 

.406 

6 

17 

50 

1 

21 

12 

el 

E2.3A 

83 

33 

09 

168 

12 

43 

146. 

. 290 

-700, 

. 982 

6 

18 

42 

1 

19 

02 

g2 


83 

33 

05 

168 

13 

14 

146. 

.211 

-701 , 
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6 

18 

4 6 

1 

18 

59 

g3 

E2 . 5 

83 

30 

04 

168 

39 

35 

141. 

.932 

-707. 

,712 

6 

22 

19 

1 

16 

40 

e 1 

E3 

83 

22 

44 

169 

34 

34 

133. 

.057 

- 723 , 

.272 

6 

30 

43 

1 

11 

53 

g2 


83 

22 

45 

169 

34 

30 

133. 

.065 

- 723 . 

,238 

6 

30 

42 

1 

11 

53 

g3 

E4 

83 

09 

20 

171 

36 

38 

110. 

924 

- 752 . 

, 137 

6 

46 

16 

0 

59 

55 

g2 


83 

09 

20 

171 

36 

40 

110. 

917 

-752. 
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6 

46 

16 

0 

59 

55 

S3 

FI 

83 

39 

43 

166 

01 

13 

170. 

048 

-683 . 

055 

6 

09 

01 

1 

31 

52 

si 


83 

39 

44 

165 

58 

18 

170. 

620 

-682. 

881 

6 

08 

56 

1 

32 

11 

s2 

F2 

83 

31 

59 
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02 

37 

173 . 
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-697. 

045 

6 

16 

34 

1 

33 

35 

s 1 


83 

31 

55 

166 

02 

57 
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-697. 
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6 

16 

38 

1 

33 

34 

s2 

G1 

•83 

41 

38 

165 

06 

18 

180. 

024 

-676 . 

816 

6 

05 

39 

1 

37 

16 

si 

G2 

83 

35 

38 

164 

23 

28 

191. 

437 

-685. 

239 

6 

10 

12 

1 

43 

25 

g3 

G3 

83 

25 

42 

162 

43 

09 

216 . 

822 

-696. 

957 

6 

16 . 

30 

1 

5 7 

08 

g2 


83 

25 

33 

162 

43 

45 

216. 

783 

-697 . 

260 

6 

16 

40 

1 

57 

06 

g3 

G4 

83 

16 

40 

161 

26 

55 

237. 

555 

- 707. 

867 

6 

22 

23 

2 

08 

19 

g2 . 


83 

16 

32 

161 

27 

17 

237. 

558 

- 708 . 

127 

6 

22 

31 

2 

08 

20 

g 3 

HI 

83 

47 

02 

164 

48 

42 

180. 

862 

-666. 

220 

5 

59 

56 

1 

37 

43 

si, 


83 

46 

58 

164 

48 

44 

180. 

888 

- 666 . 

342 * 

6 

00 

00 

1 

37 

44 

s2 

H2 

83 

46 

09 

163 

40 

13 

194. 

558 

-664. 
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5 

58 

46 

1 

4 5 

0 7 

T > 

e, ^ 

J1 

83 

35 

26 

171 

35 

33 

104. 

081 

-704. 

199 

6 

20 

26 

0 

56 

14 

P '? 


83 

35 

19 

171 

37 

00 

103. 

816 

-704. 

457 

6 

20 

34 

0 

56 

05 

h 

J 2 

83 

19 

04 

173 

04 

26 

89 . 

503 

- 736. 

795 

6 

38 

00 

0 

48 

21 

g2 


83 

18 

56 
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05 

32 

89. 
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- 73 7. 
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6 

38 

04 
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14 

, ? 
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83 

07 

22 
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55 

29 

67 . 
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- 760. 
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6 

51 

01 

0 
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30 

^2 


83 

07 

12 

174 

56 

21 

67. 

415 

-761. 
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6 

51 

11 

0 
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25 

g 3 
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83 

10 

24 
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09 

09 

155. 

683 

- 742 . 

140 

6 
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52 
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168 

09 

44 
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595 

K2 

82 
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35 
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57 

59 

136. 
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29 

169 

58 

35 

136. 
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K3 
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171 

09 

53 

122. 
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10 
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21 

101. 

022 
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04 
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59 

57 

94. 
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84 

10 

28 
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18 

43 

280. 

431 


84 

10 

28 
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21 

24 

279. 
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84 

18 

16 

154 

30 

54 

272. 

116 

A19 

84 

24 

27 

154 

42 

21 

265. 

316 


84 

24 

26 
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45 

03 

264. 
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84 

02 

42 
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07 

06 

288. 
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83 

57 

00 
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33 

32 

299. 
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58 
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35 

42 

298. 
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83 
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10 

38 
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42 
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06 

35 

290 
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12 
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E19 
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09 
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51 

07 
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53 

31 
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03 
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.813 

H2 

83 

53 

49 

150 

25 

25 

334 

. 528 


83 

53 

50 

150 

25 

29 

334 

.502 

MO 

84 

17 

46 

158 

10 

58 

235 

.383 


84 

17 

42 

158 

13 

27 

235 

.004 

Ml 

84 

22 

05 

156 

53 

27 

245 

.440 

M2 

84 

25 

40 

158 

07 

53 

230 

.456 

M3 

84 

21 

12 

159 

26 

11 

220 
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03 

g 3 

-771.946 

6 

56 

56 

1 

13 

46 

g2 

-772.152 

6 

57 

03 

1 

13 

43 

g 3 

-787.748 

7 

05 

28 

1 

06 

Q8 

g2 

-787.926 

7 

05 

33 

1 

06 

03 

g 3 

-759.397 

6 

50 

11 

0 

54 

34 

g2 

-759.429 

6 

50 

12 

0 

54 

34 

g 3 

-771.674 

6 

56 

48 

0 

51 

11 

g2 

-771.641 

6 

56 

47 

0 

51 

11 

g 3 


-582.938 

5 

14 

59 

2 

31 

31 

g2 
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2 

31 

16 
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-570.887 

5 

08 
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2 

27 

02 

el 

-561.428 

5 

03 

23 

2 

23 

21 

g 2 

-561.664 
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03 

30 

2 

23 

08 

g 3 

-594.954 

5 

21 

28 

2 

35 

58 

e 2 

-601.571 

5 

25 

02 

2 

41 

38 

g 2 

-601.815 

5 

25 

10 

2 

41 

27 

g 3 

-607.811 

5 

28 

23 

2 

46 

02 

g2 

-607.790 

5 

28 

23 

2 

46 

02 

g 3 

-572.943 

5 

09 

35 

. 2 

36 

60 

el 

-588.749 

5 

18 

07 

2 

24 

10 

el 

-597.219 

5 

22 

42 

2 

17 

58 

g2 

-597.454 

5 

22 

49 

2 

17 

45 

g 3 

-582.632 

5 

14 

48 

2 

46 

19 

g2 

-582.860 

5 

14 

56 

2 

46 

04 

g 3 

-579.350 

5 

13 

02 

2 

56 

51 

g2 

-579.533 

5 

13 

08 

2 

56 

35 

g 3 

-589.442 

' 5 

18 

28 

3 

00 

44 

g2 

-589.422 

5 

18 

28 

3 

00 

44 

g 3 

-587.986 

5 

17 

43 

2 

07 

11 

g 2 

-588.271 

5 

17 

52 

2 

06 

59 

g 3 

-575.170 

5 

10 

48 

2 

12 

38 

e3 

-574.185 

5 

10 

16 

2 

04 

32 

g 3 

-587 .036 

5 

17 

13 

1 

59 

00 

g 3 

-600.678 

5 

24 

34 

2 

01 

24 

s2 

-601.663 

5 

25 

06 

2 

09 

37 

s2 

-591.024 

5 

19 

21 

2 

14 

60 

s2 

-650.186 

5 

51 

17 

1 

54 

35 

g 3 

-669.187 

6 

01 

31 

2 

08 

55 

g 3 
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5 

43 

45 

2 

25 

01 
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5 

32 

39 
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g 3 
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APPENDIX 2 : STATIONS OCCUPIED ONCE WITH GEOCEIVERS BY NASA 
FIELD PARTIES (AS OF JANUARY 1986) 



APPENDIX 3: 


DETAILS OF MULTI-LEG ROSETTE TECHNIQUE 


Four multi-leg rosettes were planted between 1st December and 6th 
December, 1985 by two independent two-man survey teams. At each 
site, four central stakes (labeled 01 to 04) were planted into 
the snow using post levels to maintain verticality. These four 
central stakes were positioned to form the corners of a square 
with 150-m sides. Ten outlying stakes were positioned in an array 
surrounding the central square. Their distance from the central 
square was approximately 1500 m, and they were separated from 
each other by an angle of 36 degrees about 01. 

Distances were measured between each of the four central stakes 
and the ten outlying stakes using a Geodimeter 112 EDM with a 
nominal accuracy of 0.01 m over 1500 m; the times of each 
measurement were recorded to the nearest minute. Rounds of 
angles between a reference outlying stake, ' A ' and the other 9 
outlying stakes were measured using a Wild Tt 2 theodolite to a 
nominal accuracy of 3 seconds of arc, mounted on a tripod 
optically plumbed over the center of each of the four central 
stakes. The commencement time of each round was recorded. The 
true azimuth of the 'A' stake from each of the four central 
stakes was measured using observations onto the sun. The Magnavox 
MX1502 geoceiver provided both the time and the longitude used 
for the calculation of azimuth. The estimated error is 15 seconds 
of arc. This is greater than the error in observing angles 
between the outlying stakes, but the azimuth is only used to 
determine the relative vorticity (rigid body rotation) . The 
geoceiver was also used to determine the geodetic position of one 
of the central stakes to an estimated error of 20 meters, (see 
Section 1) . 

Resurvey of three of the four multi-leg rosettes was accomplished 
using the same measurement techniques by one two-man survey party 
between 16 December and 19 December, 1985. After the resurvey, 
three outlying stakes and one central stake were extended so that 
the station could be resurveyed in the future as a 3-leg. rosette . 



Data Organization 

The observed relative position of the Nth outlying stake with 
respect to each of the four central stakes is described in polar 
coordinates ( R n . 0 n . tR n . *<-» n ) where Rn denotes the radial 

distance, 6 n denotes the true azimuth (determined by adding the 
azimuth of the 'A' stake to the angular separation between the 
'A' stake and the Nth stake) , and t Rn and t Qn denote the times at 
which R n and 0 n are observed. 

The four sets of ten relative positions were combined to comprise 
one set of 40 relative positions. The effect of lumping the data 
together in this manner is to average out any strain rate 
gradients within the 150-m box formed by the four central stakes. 
This averaging is not considered detrimental because our primary 
objective is to determine the average strain rate within the 
whole rosette and not the strain rate gradient. 

Another modification made to the data was that the initial 
observations of the 'A' stake's azimuth were substituted for the 
final observations of the azimuth made during the resurvey. This 
was done to avoid introducing the azimuth uncertainty into the 
determination of stake displacements. Effectively, the azimuth 
data is used to determine the average positions of the stakes 
with respect to each other, but not to determine their relative 
displacements. Comparisons were made between strain rate 
reductions in which this substitution was and was not done, and 
no significant differences (they were less than 0.01 times the 
standard deviation of the strain rate components) were noticed. 
The relative vorticity determined as a result of this 
substitution must be corrected by adding the relative vorticity 
associated with the 'rigid rotation' determined by the rate of 
change of the true azimuth. We describe the uncertainty 
associated with the determination of the relative vorticity 
separately, below. 
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Once two sets of 40 stake positions representing the initial and 
final survey are established, average positions R n , O n and 
relative velocities R n , 0 n are computed: 


R n = y 2 (R' n + R') 


( 1 ) 


©n - 


+ 0 n) 


( 2 ) 


R„ = 


Ri - Rr 


tp - t 


(3; 


0 n = 


0 n - 0 n 

t f — t 1 


(4) 


superscripts i and f denote values determined during the 
initial and final surveys, respectively. 


These data are next converted to average positions x n , y n and 
relative velocities x n , y n in a coordinate system having an X- 
axis aligned with true North and a Y-axis aligned with true West, 


Xn = R n cos 0 n 


( 5 ) 


y n = - R n sin 0 n 
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(7) 


x n = R n cos 0 n - R n sin < ) n (-*»„ 


y n = - R n sin 0 n - R n cos B n 0 n (3) 


Determination of Strain Rates 

The determination of strain rates from a 40-leg, over-determined 
stake velocity data set requires a least-squares procedure known 
as singular value decomposition. Here we outline the details of 
this calculation. 

The strain rates are determined from the gradients in u and v 
where u and v are components of the velocity in the x and y 
directions, respectively (following Jaeger and Cook, 1976). The 
first step is to determine these velocity gradients by solving 
for S, the least-squares estimate of the velocity gradient column 
vector, in the following equation. (Note that here we write the 
equations in terms of a 40-leg rosette; the equations are easily 
modified to address an N-leg rosette where N is arbitrary.) 


AxS=d.-z (9) 


39 



where , 



and z represents the measurement error in the observed values of 
d. We assume that the covariance of z is given by: 


i = 1 ,40 C0V(Zj,Zj) 


4 COS 2 9, of, Rfsin 2 6, | 
( 4 , - t ^) 2 + < - t ;,.) 2 ) 


( 10 ) 
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1 = 41,80 cov(Z|,Zj) = 


2 . 2 a 

sin Bj 
(t' Ri - t' Ri ? 


Rfcos 2 Hi 


f >li 


where % is the Kroneker delta, CT R is the standard deviation 
estimate of the measurments of R n in meters, and is the 
standard deviation estimate of the measurement of 0 n in radians. 
This assumes covariance is consistent with the assumption that 
each of the 40 measurements of the stake velocities are 
statistically independent. However, the measurements are not 
completely independent because certain types of error associated 
with disturbing outlying stakes can affect four of the data 
points, one each associated with the observation from each 
central stake. We assume that these statistical dependencies 
produce only small off-diagonal terms in the covariance matrix of 
z compared to the diagonal terms. As stated previously, we assume 
that ff R =0.01 m and <j« =3 sec in radians. 

To correct for some components of z being larger than others y 
virtue of the inequality of a R and R n a„ , we weight the data 
matrix d and the position matrix A, 


dr = dj 


cov(z,, Zj) 


Vz 


A ?= A i 


cov(z„ z,) 


(in 


( 12 ) 


and, 
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z, 


(13) 


z 


w 

i 


cov (z,. z,) 


where the superscript w's denote weighted elements, and 


2 

a - 


80 


80 


cov (z„ z,) 

1 = 1 


(14) 


Equation (9) now becomes, 


A w x S = d w + z w 


(15) 


where the covariance matrix of z w is now a 2 on the diagonal. 

To invert A w for determining , the four eigenvalues A k and the 
eigenvectors r k of [A W ] T [A W ] are computed: 


= 


a i + 7i ( («1 + lY 


+ i : («i7i -/i?) 


16) 


+ 7i 


«1 + 7i 


4 («i7i~4 


(17) 


9 2 
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with 


and. 


40 r a 2 


40 r 

' S [; 


i = 1 


cov (Zj, Zj)J 


x 2 


(25) 


40 




I m 1 


.cov (Zj, Zj) J 


x, y, 


(26) 




40 r a 2 


.COV (Zj, Zj) J 


yf 


i = i 


(27) 


40 r ^2 -| 

2 Lcov (z i + 40 , z t + 40 )J 

j =r 1 


X? 


(28) 


4° r 

■St 


cov (z 1 + 40 , z 1 + 40 )J 


x i y, 


(29) 


40 r ^ i 

7 2 = \ 1 1 y 2 

Zj Lcov (z i + 40 , z 1 + 40 )J 1 

i = 1 


(30) 
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The least-squares estimate of S is found by constructing the 
inverse of A, 


[AT 1 = R • • Q t 


( 31 ) 


where R, r and Q are matrices defined in MacAyeal (1985) . In 
component notation, 


S; = 


- S 


k 1 

r. — 


40 


y] ( x > r i + y> r 2 ) 


k = 1 




COV (2,, Z f ) 


d. 


80 


+ 


Z <*-« r = + 

f = 41 


(32) 


where the superscript k refers to the eigenvector number. 

The principal axes el and e2 of the horizontal strain rate tensor 
are determined from S by using the relations (Jaeger and Cook, 
1976) 


0 


V2 tan 1 


/ S 2 + S 3 \ 
V S 1 + S 4 J 


(33) 
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e 1 = S, cos 2 f) + V 2 (S 2 + S 3 ) sin 20 + S 4 sin 2 0 


(34) 


and 


e 2 = St sin 2 0 - 1/2 (S 2 + S 3 ) sin 20 + S 4 cos 2 0 (3 5) 


where 0 is the counterclockwise angle (in radians) between the x- 
axis, or north, and the el-axis. 

Computation of Error 

Perhaps the most critical element of our analysis is the 
computation of expected uncertainty for the derived values of 0 , 
e , e ,and e =-e -e . This is critical for two reasons. 

1 2 ZZ 1 2 

Firstly, the short time period over which the 40-leg rosette was 
allowed to deform means that the observed stake displacements may 
be close to the limits of detectability, and secondly, the value 
of one of the strain rate principal components may be several 
orders of magnitude less than the maximum component and its value 
may be statistically insignificant. 

The expected uncertainty, or estimated standard deviation of the 
components of S , are determined formally through the singular 
value decomposition procedure (MacAyeal, 1985). 


a 


S, 


cov(Sj - Sj\ S, 



(36) 


9 6 



where S* is composed of the 7 true ( value of the velocity 
gradients, and S is, as stated before, the least-squares 
estimate of S'. Observe that the components of <1 depend only on 
the covariance of z w , the eigenvectors and the eigenvalues of L rt ] [A ], 
the latter two of which depend only on stake positioning. 

In component notation, 


ff s ( = 



(37) 


where o 2 is given by Equation (14), and subscripts i refer to 
eigenvector number. 


The uncertainty of 0, e,, e 2 are computed from the ct s .’s using the 
relationship for functions of random variables (Boas, 1983). 



(38) 


= | cos 4 9 + Va sin 2 2 0(a| 2 + a s^ + s ' n4 ® ff S 4 

+. ((S 4 - sin 2 0 + (S 2 + S 3 cos 2 0) 2 )a 2 ,,| 2 


and , 
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"e 2 = l sir|4 B lT s, + 1/4 sin 2 2 0(cr| 2 + 4 3 ) + COS 4 B a| 4 

+ ((S, - S 4 ) Sin 2 0- (S 2 + S 3 )cos 2 b) 2 <j 2 


Data reduction programs 

The above technique for reducing 40-leg rosette data can be used 
to reduce any rosette design (including the standard 3-leg 
rosette) , or other strain figures which yield data in the form 
xj. , "yi- ' x i > an< ^ Yi • A reliable calculator program which will 
reduce rosettes having up to 100 legs (a limit imposed by the 
size of calculator memory) has been developed and tested. This 
program was tested by two means. First, the results of reducing 
several 3-leg rosettes were compared with results using other 
methods, and second, multi-leg rosette data were synthesized 
using a known strain rate, and then reduced by the program to 
test whether the known strain rate is reproduced. This program, 
available on request, requires the following equipment: Hewlett- 
Packard 41-CX hand calculator having "Date" and "Time" functions 
and at least 3 modules for memory extension, an HP magnetic card 
reader to input the program, and an HP thermal printer to verify 
correct data input. The use of a hand calculator was chosen to 
allow data reduction in the field and this battery powered 
equipment operated well within the tented shelters used during 
the field program. 


Results 


Here, we present information comparing the SVD method with other 
methods of data reduction. First, the SVD method is compared to 
the method used in Section 2 of this report for calculating 
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strain rates from 3-leg rosette data to check the reliability of 
both these methods. Second, subsets of the 40-leg rosette data 
are used, which illustrate that the derived strain rate converges 
as the number of legs is increased. 

Table 18 presents strain rates derived by both methods from 3-leg 
strain rosettes surrounding the Crary Ice Rise. To illustrate 
the effect of weighting, Table 18 also includes results obtained 
when the data is not weighted in the SVD method. It is seen that 
data weighting may be considered unnecessary in the analysis of 
3-leg rosette data because errors in the distance and angle 
measurements are insignificant compared to actual changes in 
those quantities over the survey period. For all 14 stations, the 
principal horizontal strain rate components derived by the two 
methods differed by less than 10 percent and fall within the 
computed uncertainty limits obtained by both methods. This 
comparison serves as a useful check on both our data-analysis 
techniques . 

To illustrate convergence of the SVD technique in producing an 
accurate estimate of the horizontal strain rates from the 40-leg 
rosette N3, twelve independent 3-leg rosettes and four 
independent 10-leg rosettes from the 40-leg rosette data are 
synthesized by considering subsets of the stake array. The 
results are given in Table 18. There is considerable scatter in 
the results of the twelve 3-leg rosette as expected from the 
short time period before which the 3-leg rosettes were 
resurveyed, but the strain rate falls within the sample standard 
deviation of the strain rate derived from the analysis of the 40- 
leg data. The sample standard deviation of the twelve 3-leg 
rosettes results is larger than the confidence limits computed 
from 40-leg rosette data. This could have resulted from three 
factors: 1) our estimates of measurement error on R and 0 are 
too low, 2) the sample of twelve 3-leg rosettes is too small to 
produce an accurate sample standard deviation, and 3) there is 
strain rate variation within the confines of our stake array that 
is averaged out by the 40-leg analysis. All these possible 
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explanations have some validity. However, the possible 
underestimate of the measurement accuracy seriously affects the 
scientific conclusions of our study; certainly any underestimate 
is not a factor of 10, which is needed if that alone causes the 
discrepancy . 

As another illustration of the multi-leg rosette analysis, four 
10-leg rosette sub-sets were analyzed, each of the 4 central 
stakes being the center of a separate stake array. There is much 
less scatter in these results (Table 18) when compared to the 
twelve 3-leg rosettes. The sample standard deviation is also 
consistent with the confidence limits of the 40.-leg rosette 
showing there is a convergence in both the derived strain rate 
and the confidence limits as the number of stakes is increased. 


INSTRUMENT CORRECTIONS 


Benchmark comparisons were made periodically between all EDM's 
used in the 1985-1986 field season. These comparisons revealed 
that 1 EDM used in the initial and final surveys of 40-leg 
stations N1 and N4 requires a -5 . 5xl0~ 5 +0 . lxlO 5 parts-per-part 
correction to the initial measured distance. 


MULTI-LEG ROSETTE UTILITY 


Multi-leg rosettes can serve a useful purpose under circumstances 
when more practical 3-leg rosettes will not provide sufficient 
accuracy over the time interval between survey and resurvey. 
Consideration must be taken of the large commitment in time 
required to deploy these rosettes. Experienced surveyors, in good 
weather, needed 1.5 days to deploy a 40-leg rosette and 1 day to 
resurvey it. The method also requires instrument dependability. 
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